Jamming of 3D Prolate Granular Materials 
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We have found that the ability of long thin rods to jam into a solid-like state in response to a local 
perturbation depends upon both the particle aspect ratio and the container size. The dynamic phase 
diagram in this parameter space reveals a broad transition region separating granular stick-slip and 
solid-like behavior. In this transition region the pile displays both solid and stick-slip behavior. We 
measure the force on a small object pulled through the pile, and find the fluctuation spectra to have 
power law tails with an exponent characteristic of the region. The exponent varies from j3 = —2 
in the stick-slip region to /? = — 1 in the solid region. These values reflect the different origins — 
granular rearrangements vs. dry friction — of the fluctuations. Finally, the packing fraction shows 
only a slight dependence on container size, but depends on aspect ratio in a manner predicted by a 
mean field theory and implies an aspect-ratio independent contact number of (c) = 5.25 ± 0.03. 



I. INTRODUCTION 

Anyone scooping nails from a bin at the hardware store 
or shoveling hay with a pitchfork has noticed the ability 
of large aspect-ratio (a = L/D ^> 1) granular materials 
to act as a solid. This state can exist at low packing frac- 
tions, raising the possibility of lightweight, rigid building 
materials. Such materials have practical applications as 
proppant in oil recovery or conducting networks in quan- 
tum dot solar cells. Here we present the first experi- 
mental investigation into the dependence of this state on 
particle aspect ratio and container size. 

The jamming of round and low aspect-ratio (< 5) ma- 
terials has attracted a great deal of recent attention. A 
canonical example of granular jamming is flow through 
hoppers, where arches at the exit can support the entire 
pile 0, 0, El H- This jammed state, while capable of 
supporting large constant forces, is not robust when sub- 
jected to sharp taps. Indeed, in one set of experiments 
4] a jet of air was sufficient to reintroduce flow. 

Jamming is also invoked to explain the erratic drag 
force on an object moving through a granular medium 
H El. The object comes to rest when it encounters a 
connected network of particles that is incapable of rear- 
ranging; only when the applied force is large enough to 
break through this network does the motion begin again. 
The particles in this network are only a small fraction 
of the total grains in the pile, those comprising the force 
chains |7j,|8j that terminate on the container walls. 

The situation is remarkably different for high aspect- 
ratio (> 10) granular materials Philipse [l(| noted 
that particles of aspect ratio larger than 35 emerged as a 
solid plug when poured. This plug maintained the shape 
of the original container even in face of external distur- 
bances. This jammed state can therefore be considered 
both global and robust, the connected network of parti- 
cles containing a significant fraction of the pile that are 
incapable of moving around one another. How this solid- 
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like state compares with the jammed state of ordinary 
granular materials has not yet been explored. 

Recent work has investigated the mean drag force ex- 
erted by ordinary granular materials on an intruder. 
Zhou et al. found [ill Il2|. for example, that mean drag 
scales linearly with pressure for both mono- and poly- 
disperse materials. Geng and Behringer^l studied the 
force fluctuations exerted by a two-dimensional packing 
of disks and found that the power spectra had a power 
law tail with exponent —2, as in three-dimensional sys- 
tems. Hill et al.0 and Stone et al.[D| 13 have all 
looked at plates pushed slowly down into beds of sand, 
with interesting behaviors as the plate approaches the 
bottom of the container and a temporal evolution in the 
granular bed with successive plungings. Finally, Brat- 
berg et al.^3 simulated the quasi-static flow of rigid, 
frictional disks pushed upward against gravity through 
a narrow pipe. They observed a transition in the flow 
when the intruder speed was large. 



II. SETUP AND PACKING FRACTION 

In our experiments, the force on a small ball pulled 
upward through a pile of prolate granular materials is 
measured. Particles are cut from acrylic rods of diam- 
eter d — 1/8" or d — 1/16" to varying lengths; aspect 
ratios vary from 4 to 48. The particles are dropped into 
cylinders with diameters ranging from D = 1" to D = 4". 
A dimensionless container diameter is formed by divid- 
ing by the particle length (D = D/L); the smallest D 
is 0.67, the largest 12. At the bottom of the pile sits a 
small (d = 0.25") metal ball connected to string that runs 
up through the pile. We have confirmed that packings 
formed around the string have the same packing fraction 
as those formed in the string's absence. The string is 
wound around the shaft of a motor that turns with con- 
stant angular speed. A force sensor measures the tension 
in the string to within 0.5% at kHz resolution. For the 
range of forces in our experiments, the string and force 
sensor act Hookian with spring constant k = 1600 N/m. 
The motion of the ball can be erratic; when the resis- 
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tance of the pile exceeds that applied by the string the 
ball will stop moving until the applied force, increasing 
as the string is wound, is again large enough to induce 
motion. 

Particles are simply dropped into the cylinder. "Fluff- 
ing" the pile with air, a technique used to prepare ordi- 
nary granular materials in a low density state, does not 
significantly decrease the packing fraction of rods [l8| . 
The packing fraction is plotted (log-log) versus aspect 
ratio in Fig. ^ with different symbols indicating different 
dimensionless container sizes. The data from different 
size containers are reasonably close to one another, ex- 
cept for a few points at the low D. At these values the 
particles tend to align with the sidewalls, introducing sig- 
nificant orientational correlation and higher packing frac- 
tions. 
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FIG. 1: Log-log plot of packing fraction <f> as a function of 
particle aspect ratio a. Data are grouped by constant dimen- 
sionless container diameter D = D/L. With the exception of 
a few points at low D (where significant particle alignment 
occurs), the data are well- fit by a mean- field theory with only 
one free parameter. 

To gain some idea of the randomness of the initial pack- 
ing we compare the packing fraction with a mean-field 
theory the Random Contact Model (RCM) [Tjj. The ba- 
sic assumption of this model is that contacts between 
particles are uncorrelated, a significant difference from 
the correlated contacts that result in ordered packing of, 
for example, disks or rods. If the contacts are uncorre- 
lated, then the packing can be thought of as a collection 
of independent pairs of rods in contact. 

A key finding of Philipse is that the average number 
of contacts (c) scales linearly with with the particle con- 
centration. To show this, Philipse begins with the frac- 
tion of orientations made inaccessible to a particle by the 
existence of another particle. This is a function of the 
center-of-mass separation of the two particles; when two 
rods are close together, they must be more nearly aligned 
to prevent overlapping. Earlier work on the packing of 
rods in two dimensions |9| calculated the functional de- 
pendence of the excluded fraction of orientations / ex (r) 



analytically and found that, while both experiment and 
simulations contained long-range correlations, the func- 
tional form was universal over several orders of magni- 
tude in aspect ratio. The average number of contacts is 
then proportional to the product of the excluded fraction 
of orientations and the local number density: 

(c) = \ J /ex(p, r)p(f)df. 

The proportionality factor of 1 /2 avoids double counting 
of contacts. 

The mean-field approximation replaces the local den- 
sity p{r) with the average pile density p; the remaining 
integral over the excluded orientations can then be iden- 
tified with the average excluded volume v eyi , defined 0] 
as the volume denied to particle j by the condition that 
it not overlap with particle i. The average contact num- 
ber then scales linearly with packing concentrations and 
the number density of a pile is 
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The excluded volume for cylinders of length L and 
thickness D is 191 
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The predicted packing fraction 4>(a) has one single free 
parameter, the average contact number (c). We (and 
others |18|') noticed that significant orientational order 
exists when the cylinder diameter D is less than 1 par- 
ticle length L. This violates the main assumption of 
the model, and so we should not expect good agreement 
in this region. We therefore attempt to fit the predic- 
tion from Eq. [3 with all data taken in cylinders where 
D/L > 1. The resulting line is shown in Fig. 2] For 
large aspect ratios Eq. |21 reduces to 6 ~ (c) /a, a scal- 
ing that was noticed previously in 2-[9j and 3-dimensions 
(l0|. The curvature in the data in Fig. ^ shows that this 
simple limit does not apply; the full expression of Eq. [21 
(solid line) agrees quite well with the data. From the fit 
we find (c) = 5.25 ± 0.03, which agrees (within reported 
error) with the data reported by Philipse . The con- 
stant contact number implies a significant screening ef- 
fect, since the contact number could in principle scale 
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with aspect ratio. A similar result for two dimensional 
piles was found in earlier simulations [jj- 

The good agreement between experiment and a model 
that assumes an isotropic distribution of particle angle 
supports the claim that particles are initially randomly 
oriented. The discrepancy at small container size is un- 
derstandable, as the orientational correlation induced by 
the narrow cylinder violates the fundamental assumption 
of the RCM. 



III. FORCE SCALING IN DIFFERENT 
REGIONS 

When the particle aspect ratio is low, the pile re- 
sponds with local rearrangements and the drag force on 
an intruding object has a random-sawtooth appearance 
(Fig-IHtop)). The force increases linearly, indicating that 
the ball is at rest, before rapidly decreasing, indicating a 
burst of motion. The ball is brought to rest and the cycle 
repeats. Throughout the experiment, the bulk pile is at 
rest. While individual particles near the ball are moving, 
there is no collective pile motion. As the aspect ratio is 
increased, however, the pile exhibits a qualitative change 
to solid-like behavior. With the exception of a few stray 
particles, all particles are lifted upward, with no observed 
relative motion between particles. The resulting force vs. 
time diagram is shown in Fig. [^bottom) . Forces are nor- 
malized by the total pile weight, and so Fig. [21 shows that 
the force required to move the pile can be many times the 
actual pile weight. This is a consequence of force chains 
that terminate on the container walls. The ball is 
pushing upward on the particles, but the lateral deflec- 
tions of the force chains translate this into a force normal 
to container walls. This normal force can be quite large, 
resulting in large frictional forces between the walls and 
the pile. 

For intermediate aspect ratios, the pile displays char- 
acteristics of both the granular and solid states. Fig- 
ure |2I middle) shows data from one such run. The pile 
first behaves much like that of smaller aspect ratio par- 
ticles, although a close examination of the force (inset 
in Fig. |2I middle)) reveals small plateaus indicating small 
amounts of solid-body behavior. Then, with the ball in 
the middle of the pile, the entire pile jams and all par- 
ticles above the ball are lifted. Corresponding to the 
plateaus in the force is the collective motion of the entire 
pile. The pile is visually observed to move as a solid body 
for a brief period of time. The breakup of the solid state 
is brought about by a collapse of the particles around the 
ball. As the pile is constantly rubbed against the side- 
wall (and, indeed, shows force fluctuations characteristic 
of dry friction) the length of time spent in the solid state 
is some indication of the pile's stability. 

A Fourier transform of the force vs. time data in the 
granular and solid regions (Fig. [3J shows that both spec- 
tra show power law tails, albeit with different exponents 
indicating a different fluctuation origin. Data from the 
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FIG. 2: Force vs. time data for a ball dragged through a pile 
of rods of aspect ratio 12 (top), 20 (middle), and 40 (bottom), 
all in a 2" tube. The low aspect ratio particles show the 
stick-slip behavior common in ordinary granular materials. 
The large aspect-ratio particles, however, act as a single solid 
body, with small fluctuations characteristic of dry friction. 
Intermediate aspect ratios show both behaviors in a single 
experiment on both large and small (inset) time scales. 



granular region decay as P(f) ~ /~ 2 , consistent with ear- 
lier experiments on granular materials [2£j, El • Coupled 
with the visual observation of no macroscopic pile mo- 
tion, this supports connecting the force serrations with 
local particle rearrangements. The forces in the solid re- 
gion, however, decay as P(f) ~ / . This is consistent 
with experimental work on dry friction |22j |. and so we 
infer that the pile is not moving steadily upward, but in 
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FIG. 3: Power spectra for force fluctuations from piles exhibit- 
ing granular (o) and solid (x) behavior. Both show power law 
tails, with the different exponents indicating a different mech- 
anism for the fluctuations. The 1// and and l// 2 decay are 
consistent with previous work on, respectively, dry friction 
and localized granular rearrangements. Both data sets taken 
in a D=2" diameter tube. 

fact sticking on the side walls. 

The exponent of the spectrum tail is a reliable indi- 
cator of the pile's behavior. Shown in Fig. 0] is the ex- 
ponent magnitude | j3 | as a function of normalized tube 
diameter D. The data are grouped in three sets corre- 
sponding to the three different behaviors: stick-slip (o), 
transition (o), and solid- like (□). Lying off the graph's 
scale are three additional stick-slip data points (each rep- 
resenting an average of several runs) at D = 5, 8, &12; all 
are characterized by an exponent of magnitude of just un- 
der 2 with relatively small uncertainties. Data exhibiting 
transitional behavior show much larger uncertainties, in- 
dicating the statistically fluctuating nature of the behav- 
ior. While all force spectra from these experiments show 
the plateaus characteristic of the transition region, the 
relative amount of time spent on a plateau (as opposed 
to stick-slip) varies wildly from run to run. Data in the 
solid region have exponents close to one with small un- 
certainties (with the exception of data at D — 1.2, which 
we suspect was very close to the transition region). 



IV. DYNAMIC PHASE DIAGRAM 

To better understand the transition from granular to 
solid-like motion, we have mapped out the dynamic phase 
diagram of the behavior as a function of the two control 
parameters — aspect ratio a and container diameter D. 
It is in fact more revealing to use the inverse container 
diameter 1/D. This is because when the particle length 
is increased but tube diameter held fixed both the aspect 
ratio and inverse diameter (<5 = L/D) increase linearly, 
and so a sequence of experiments in which particle length 



FIG. 4: Exponents of the spectra tails for experimental runs 
of varying aspect ratios plotted vs. normalized container di- 
ameter D = D/L. Filled symbols are aspect ratios larger than 
20, which display solid-like motion in containers less than 1.5 
particle lengths across. 

is increased shows up as straight lines in S— a space. This 
is shown in Fig. [5] 
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FIG. 5: Phases exhibited by granular piles as a function of 
two control parameters — the aspect ratio a and the inverse 
container diameter 5 = D . Smaller aspect ratio particles 
show the stick-slip behavior of granular materials, while larger 
aspect ratio particles act as a solid body when the container 
is small enough. 

As one would expect, when the aspect ratio is very 
small, the pile behaves in a canonically granular manner. 
Nevertheless, the signature characteristics of the transi- 
tion region — plateaus in the force data and a visual 
observation of collective motion — are seen in particles 
with aspect ratios as low as 8 when confined to cylin- 
ders whose diameter is twice the particle length. As the 
container diameter is further reduced, and the ability of 
particles to move around one another further constrained, 
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we expect solid-body motion to occur even for the small 
aspect ratios (upper left corner of the phase diagram), 
but we did not explore this region. 

We did not investigate the dependence of the transi- 
tion region on pulling speed. There are two different 
mechanisms through which the pulling speed can be- 
come significant: elastic unloading (particularly signifi- 
cant in a stiff machine) and particle rearrangements. We 
show here that we are well below the critical velocities at 
which cither of these become important. In a theoretical 
investigation of frictional disks pushed upward through 
a vertical pipe, Bratberg et al.|l7j found a transition 
that depended on the relative value of two time scales, 
tk = (x) 1 ^ 2 ( TO tne P ne weight, k the spring constant 
of the pulling mechanism) corresponding to the elastic 
unloading and t g = - the time for a particle to acceler- 
ate due to gravity to velocity v. An order-of-magnitudc 
estimation of these values for our system is m ~ 1 kg, 
v = 0.01 m/s, and k = 1600 N/m, and therefore 

t g v/g ~ 10~ 3 s ~ 

Using these numbers, we would expect a significant in- 
fluence of pushing velocity to occur at a characteristic 
velocity of about v c = \prQ ~ 24 cm/s, twenty times 
our actual pulling speed. Therefore, our pulling "spring" 
is soft enough, and our velocities slow enough, so that 
the drop in force is not significant. 

If the particle rearrangements are governed solely by 
gravitational forces, then the ratio of the time it takes a 
particle to fall under gravity a distance equal to its length 

L t r = \ — to the time it takes the intruder to move that 

distance ti = L/v is important. The ratio of these times, 
using our approximate parameters, is then 

^ = ^«50, 

t r V 

leading to a critical velocity of v c = \fg~L « 50 cm/s, 
about 50 times faster than the current pulling speed. 
Thus, the intruder moves a negligible distance while the 
grains themselves are re-orienting. The above analyses 
imply that we are in the quasi-static region of behavior 
and will not see significant changes in the behavior un- 
less the pulling speed is increased by at least an order of 
magnitude. 

Philipse observed that particles with aspect ratios 
larger than 35 did not flow when poured from a bucket, 
but rather emerged as a single, solid plug. The expla- 
nation given was a geometric entanglement of the rods. 



Interestingly, however, we do not observe the solid behav- 
ior in particles with an aspect ratio of 48 when poured 
into the large cylinder. It may be that, in response to 
the localized disturbance of a small intruder, the pile 
can make small rearrangements necessary to allow the 
intruder to pass through while still maintaining an over- 
all rigid structure. An important next step in this work 
will be to investigate the effect of intruder size on phase 
behavior. For small intruders, we imagine we are prob- 
ing details of the spaces (voids) in the pile while larger 
intruders probe the pile's rigidity. This would imply a 
critical length scale related to particle length and width, 
consistent with our observation that the particle aspect 
ratio alone does not completely determine the pile's be- 
havior. 



V. CONCLUSIONS 

We have observed three qualitatively different types of 
behavior in large aspect-ratio granular materials in re- 
sponse to a local disturbance: canonically granular stick- 
slip, solid-body like motion, and a transition region that 
is a combination of the two. The phase space of this 
behavior has been mapped out as a function of two con- 
trol parameters, the particle aspect ratio and container 
diameter. We have found that even low aspect-ratio par- 
ticles can exhibit temporary solid-body motion that is 
characteristic of transitional behavior. Surprisingly, the 
larger aspect ratio particles do not behave as solid bod- 
ies when the container is large enough. The granular 
and solid states both show force fluctuations with spec- 
tra that have power-law tails, although the exponents are 
characteristically different. 

We have also investigated the packing fraction as a 
function of these two control parameters. While there is 
a slight dependence on container size, we find a distinct 
dependence on aspect ratio that agrees quite well with 
the mean-field Random Contact Model. This is perhaps 
surprising, as we observe some orientational ordering that 
violates a main assumption of this model. 
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